NMR LOGGING USING TIME-DOMAIN AVERAGING 



FIELD OF THE INVENTION 

5 The present invention concerns logging applications and more specifically relates to 

a method for data acquisition and processing, which uses time-domain averaging to increase 
the signal to noise ratio (SNR) or to sharpen the spatial resolution of the obtained logging 
data. In an exemplary application, the method is used in nuclear magnetic resonance 
(TSfMR) logging to increase the SNR of single event echo-train for high-speed NMR logging 

10 and/or to sharpen the spatial resolution of the received data. 

BACKGROUND OF THE INVENTION 

In oil and gas exploration it is desirable to understand the structure and properties of 
the geological formation surrounding a borehole, in order to determine if the formation 

15 contains hydrocarbon resources (oil and/or gas), to estimate the amount and producibility of 
hydrocarbon contained in the formation, and to evaluate the best options for completing the 
well in production. A significant aid in this evaluation is the use of wireline logging and/or 
logging-while-drilling (LWD) measurements of the formation surrounding the borehole 
(referred to collectively as "logs" or "log measurements"). Typically, one or more logging 

20 tools are lowered into the borehole and the tool readings or measurement logs are recorded 
as the tools traverse the borehole. These measurement logs are used to infer the desired 
formation properties. 

In recent years nuclear magnetic resonance (NMR) logging has become very 
important for purposes of formation evaluation and is one of the preferred methods for 

25 determining formation parameters. Improvements in the NMR logging tools, as well as 
advances in data analysis and interpretation allow log analysts to generate detailed reservoir 
description reports, including clay-bound and capillary-bound related porosity, estimates of 
the amounts of bound and free fluids, fluid types (i.e., oil, gas and water), permeability and 
other properties of interest. 

30 NMR tools used in practical applications include, for example, the centralized 

MRIL® tool made by NUMAR Corporation, a Halliburton company, and the sidewall CMR 
tool made by Schlumberger. The MRIL® tool is described, for example, in U.S. Pat. 
4,710,713 to Taicher et al. and in various other publications including: "Spin Echo 
Magnetic Resonance Logging: Porosity and Free Fluid Index Determination," by Miller, 

35 Paltiel, Gillen, Granot and Bouton, SPE 20561, 65th Annual Technical Conference of the 
SPE, New Orleans, LA, Sept. 23-26, 1990; "Improved Log Quality With a Dual-Frequency 
Pulsed NMR Tool," by Chandler, Drack, Miller and Prammer, SPE 28365, 69th Annual 
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Technical Conference of the SPE, New Orleans, LA, Sept. 25-28, 1994. Certain details of 
the structure and the use of the MRIL® tool, as well as the interpretation of various 
measurement parameters are also discussed in U.S. patents 4,717,876; 4,717,877; 
4,717,878; 5,212,447; 5,280,243; 5,309,098; 5,412,320; 5,517,115, 5,557,200; 5,696,448; 

5 5,936,405; 6,005,389; 6,023,164; 6,051,973; 6,107,796 and 6,111,408. The structure and 
operation of the Schlumberger CMR tool is described, for example, in U.S. Pat. Nos. 
4,939,648; 5,055,787 and 5,055,788 and further in "Novel NMR Apparatus for 
Investigating an External Sample," by Kleinberg, Sezginer and Griffin, J. Magn. Reson. 97, 
466-485, 1992; and "An Improved NMR Tool Design for Faster Logging," D. McKeon et 

10 al., SPWLA 40 th Annual Logging Symposium, May- June 1999. The content of the above 
patents is hereby expressly incorporated by reference for all purposes, and all non-patent 
references are incorporated by reference for background. 

NMR tools of the type discussed above generally measure the time for hydrogen 
nuclei present in the earth formation to realign their spin axes, and consequently their bulk 

1 5 magnetization, either with an externally applied magnetic field, or perpendicularly to the 
magnetic field, after momentary reorientation due to the application of specific radio 
frequency (RF) pulses. The externally applied magnetic field is typically provided by a 
magnet disposed in the tool. The spin axes of the hydrogen nuclei in the earth formation 
are, in the aggregate, caused to be aligned with the magnetic field induced in the earth 

20 formation by the magnet. The NMR tool includes an antenna positioned near the magnet 
and shaped so that a pulse of radio frequency (RF) power conducted through the antenna 
induces a magnetic field in the earth formation orthogonal to the field induced by the 
magnet. The RF pulse has a duration predetermined so that the spin axes of the hydrogen 
nuclei generally align themselves perpendicular both to the orthogonal magnetic field 

25 induced by the RF pulse and to the externally applied magnetic field. After the pulse ends, 
the nuclear magnetic moment of the hydrogen nuclei gradually relax, i.e., return to their 
alignment with the externally applied magnetic field; at the same time an antenna, which is 
typically the same as the one used by the initial pulse, is electrically connected to a receiver, 
which detects and measures voltages induced in the antenna by precessional rotation of the 

30 spin axes of the hydrogen nuclei. 

An actual NMR measurement involves a plurality of pulses grouped into pulse 
sequences, most frequently of the type known in the art as Carr-Purcell-Meiboom-Gill 
(CMPG) pulsed spin echo sequences. As known in the art, each CPMG sequence consists 
of a 90-degree (i.e., tz/2) pulse followed by a large number of 1 80-degree (i.e., it) pulses. 

35 The 90-degree pulse rotates the proton spins into the transverse plane and the 1 80-degree 
pulses generate a sequence of spin echoes by refocusing the transverse magnetization after 
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each spin echo. A wait time interval (T w ) of approximately 0.5-10 sec is used first to allow 
for polarization of the formation by the tool's static magnetic field. Then, a CPMG 
pulse-echo train is executed, consisting of an excitation pulse and an alternating sequence of 
refocusing pulses. Following each pair of excitation pulse and a refocusing pulse, an 

5 acquisition window is applied next. Complex data points from such a pair of echo trains are 
co-added on an echo-by-echo basis to remove certain artifacts and to enhance the NMR 
signal, as known in the art. More pairs may be added to enhance the signal-to-noise ratio. 
The echo train, consisting of a superposition of exponentially decaying signals is then 
submitted to a processor, which calculates the underlying decay modes of the received 

10 NMR echo signal. 

One of the key information sought from the NMR logging data is the T 2 decay time 
distribution, which is related to the pore size distribution of the rock formation. From the T 2 
decay time distribution one can estimate, for example, the capillary bound water volume. 
The T 2 distribution is obtained by inversion of the echo train data, using well-known 

15 methods as described, for example, in U.S. Pat. No. 5,517,1 15 to the assignee of the present 
application. 

It is also well-known in the art that inversion of NMR echo trains is greatly affected 
by noise. Fig. 1 illustrates a single-event echo train. In order to reduce the effect of noise 
and increase the signal-to-noise ratio (SNR), stacking of multiple echo train data is 

20 commonly used, by assuming that noise is of random nature. Thus, for example, the 

average echo train from 10 events is affected with much less noise than the echo train from 
a single event. (In general, if N is the number of stacked signals, the SNR of the averaged 
signal improves as approximately sqrt(N)). Various methods exist to perform stacking 
such as, for example, using two pulse sequences of opposite phase that cancel electronic 

25 offsets and 180-degree ringing. Such a pair of pulse sequences is called a phase-alternated 
pair (PAP). 

In actual logging, the tool is moving at a certain speed, and echo train from each 
event is not collected at the same logging location or depth. When the logging speed is 
slow, the tool moves only a short distance. Hence one may assume that the earth formation 

30 is of nearly constant property over the short distance, and stacking of multiple events from 
different depths may be acceptable. However, when the logging speed is relatively fast, the 
tool moves over longer distances and the formation properties may change significantly 
over the distance. As a result, stacking of echo train data from multiple events, as generally 
used in the prior art, reduces apparent spatial resolution of measuremenets beyond the 

35 intrinsic resolution of the tool. 
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In particular, in high-signal formations (e.g., high porosity, oil-or water-filled 
rocks), where the logging speeds can be comparatively fast, the effective vertical resolution 
of the NMR log becomes dominated by the tool movement during a single PAP. For 
example, with a recovery time of 1 0 seconds between echo-trains in a PAP, with P APs 

5 acquired at all possible frequencies, the elapsed time between the first echo-train in the first- 
frequency PAP and the second echo-train for the last-frequency PAP, is close to 20 seconds. 
At a logging speed of 900 ft/hr (15 ft/min), the MRIL tool will move approximately 5 feet 
during this measurement: when combined with the inherent vertical resolution of the 
antenna (which is approximately 2 feet), the effective vertical resolution becomes roughly 7 

10 feet. 

Enhancing the resolution of the logs is a significant problem, because subsurface 
formations are generally heterogeneous, so that porosity, saturation and lithology vary with 
position. A common example of heterogeneity is the presence in the formation of 
geological layers, or beds. Because logging tools have a nonzero volume of investigation, 

1 5 more than one layer may lie within the volume of investigation of a tool. In such cases, the 
petrophysical evaluation of one layer may be distorted by the presence of another layer 
falling within the larger volume of investigation of the tool. The above phenomenon leads 
to a specific problem in the analysis of subsurface formations that include one or more 
underground layers, especially when the layers are thin compared with the vertical 

20 resolution of the measuring tool. Such layers have become subject to significant 
commercial interest because of their production potential. Any knowledge about the 
composition and properties of such layered formations that helps better estimate their 
production potential has thus become increasingly valuable. It is clear therefore that any 
mechanism that for a given SNR supported by the formation can increase the vertical 

25 resolution of the tool without decreasing the logging speed is highly desirable. 

SUMMARY OF THE INVENTION 

In a preferred embodiment the present invention uses time-domain average of single 
event data to increase the SNR, but not stacking of data from multiple events, and hence, 

30 achieves higher spatial resolution. In accordance with the present invention, the SNR can 
be increased also by application of both the time-domain-average and stacking techniques. 
In such case, high spatial resolution can be achieved using a lesser number of stacking 
events. Various combinations of techniques known in the art with the teachings of the 
present invention are possible in this regard and will be appreciated by a person of skill in 

35 the art. 
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In particular, according to this invention a method for measuring an indication of 
attributes of materials containing a fluid state is disclosed, the method comprising the steps 
of: providing a time-domain signal indicative of attributes of said materials in a single 
measurement; constructing a time-domain averaged data train from said signal, the 
5 averaging being performed over one or more time intervals and computing an indication 
of attributes of said materials from the time-domain averaged data train. In different 
embodiments, the time intervals A { can be selected constant or different from each other. 
In a specific preferred embodiment of the method the time-domain signal is an NMR echo 
train. 

10 In accordance with another aspect of the invention, a method for increasing the 

spatial resolution of NMR logging measurements is disclosed, comprising the steps of: 
providing an NMR echo-train indicative of attributes of materials of interest; and 
constructing a time-domain averaged data train from said NMR echo train, the averaging 
being performed over one or more time intervals It should be understood that the 

15 method of the present invention is applicable to processing of data from virtually any type 
of logging device and thus is not limited to NMR logging applications. 

BRIEF DESCRIPTION OF THE DRAWINGS 

20 Fig. 1 is an illustration of a single event echo train from a NMR log; 

Fig. 2 illustrates a model T 2 distribution and inverted T 2 distribution obtained from a 
single event echo train; 

Fig. 3 illustrates multi-event echo stacking; 
Fig. 4 is an illustration of a stacked echo train from ten events; 
25 Fig. 5 illustrates an inverted T 2 distribution from stacked echo-train data; 

Fig. 6 illustrates time-domain averaging used in accordance with the present 
invention; 

Fig. 7 illustrates a time-domain averaged echo-train data; and 

Fig. 8 illustrates an inverted T 2 distribution from time-domain averaged echo-train 



30 data. 



35 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The method in accordance with the present invention will be described next in the 
context of NMR logging and will use terminology generally accepted in this field. It should 
be understood, however, that the method is equally applicable to processing of data from 
any number of other logging applications, and more generally to processing of data that can 
be modeled as an exponential function of time. 

With reference to the NMR logging applications discussed above, in accordance 
with the present invention, the T 2 decay time distribution can be obtained by inverting the 
time-domain echo train data, assuming the following physical model: 



10 



5(0 = X ^(^) ex P(-' fT 2 )+ Noise 



20 



where <j) (T 2 ) is the porosity (or population) of the pores corresponding to the exponential 

15 decay time T 2 . One such example of echo-train is plotted in Fig. 1, where the echo train is 
generated from the T 2 distribution, $ (T 2 ), which is illustrated in Fig. 2. Because of noise, 

the inverted T 2 distribution obtained from this single event echo-train, which is also shown 
in Fig. 2, is different from the model input. 

To reduce the effect of random noise, it is commonly practiced to stack echo-trains 
from a mult^te^iuiaberof events. Fig. 3 illustrates the basic idea behind stacking. The top 
portion of the figure illustmtB^aglo^gvent echo train obtained from a moving tool having a 
given intrinsic vertical resolution deteiminedpmtma^by the dimensions of the tool. The 
bottom portion of the figure shows the effect of multi-evenT^ekiiig, which clearly reduces 
the noise considerably, but also acts to change the apparent vertical resolttttea^of the to ol. 

Fig. 4 shows stacked, i.e., averaged, echo trains obtained from a combination of 10 
events. Significant reduction of random noise is observed, as indicated by the smooth ten- 
event average line. The inverted T 2 distribution from the 10 event stacked echo train is 
closer to the model input, as illustrated in Fig. 5. 

As indicated in Fig. 3, stacking of multiple event data causes worsening of the 
spatial resolution of NMR logging, as the tool is moving in usual logging measurements. 
As noted above, the tool moves a considerable distance at high logging speeds. It should be 
apparent that the formation property may change significantly over the distance. While 
single-event echo-train data provides information of the formation properties for an interval 
determined by the intrinsic spatial resolution of the tool, the stacked echo train data 
provides formation property average over the distance to cover such multiple events. As a 
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result, the spatial resolution of the measurement is greatly reduced beyond the intrinsic 
resolution of the tool. 

In accordance with the present invention, the signal-to-nose ratio (SNR) can be 
increased without stacking echo train data from multiple events (or by using reduced 
stacking) and thus without worsening the spatial resolution at high logging speeds, but by 
taking time-domain-average(s) of the input signal. The principal advantage of this approach 
is therefore that the spatial resolution of the tool over the coarse of a logging experiment is 
no-Lworsened at the price of the resolution in the time domain. 



-domain approach to enhancing vertical resolution of data logs in 
10 accordance with the^presgnt invention is illustrated in Fig. 6. In particular, in accordance 
with the present invention theT^k^ing approach is used. First, construct a time domain 
averaged echo train as follows: for time vfafee^ao^yerage of echos over a time interval A, 
as defined by the following expression: 



t + A 



15 



S A (0 = jdt' S(t' ) / A 



It is important to note that the noise in the averaged echo is still random. To ensure 
that, the averaged echo separated by time interval of A is considered. Namely, the averaged 

2^ echo train at t= t 0 , t 0 + A, t 0 + 2 A, .... , t 0 + NA is used to obtain T 2 distribution. Fig. 7 shows 
such a time-domain averaged echo train averaged over 10 sampling points. In this example, 
the single echo train has 300 data points at t/T E = 1, .... , 300 (with T E = 1.2 ms). As A = 12 
ms (10 sampling points), the averaged echo has 30 (=300/10) data points. The inverted T 2 
distribution obtained from this time-domain-averaged echo is closer to the model input than 

^ the one from a single echo, as illustrated in Fig. 8. Note that the following equation is used 
to estimate T 2 distribution: 

(0 = X 0( T i ) ex P(" * I T 2 )(1 - exp(- A / T 2 )) + Noise 

It should be apparent that in different applications the value of the parameter A can 
30 be changed to conform with the underlying physics of the signal. Furthermore, in 

alternative embodiments of the invention, instead of a single time period A, one can use two 
or more periods A { the values of which can generally be different. As in the preceding case, 
the values for Aj can be selected dependent on the practical application and in particular on 
the need to increase the SNR of the received signal. 

35 
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In accordance with another aspect of the present invention, stacking of multiple- 
event echo-train data and time-domain averaging can be used in combination. For instance, 
instead of taking 10-event stacking or time-domain averaging over 10 sampling points, one 
can use 5 -event stacking and then apply time-domain averaging to increase the SNR of the 

5 signal, without worsening the vertical resolution of the log as much as prior art methods 
would, which only use stacking. 

It should be noted that while the method of the present invention has been described 
in the context of NMR logging, it is applicable to a multitude of practical applications 
which need not be limited to NMR logging. Virtually all types of logging measurements 

10 with relatively low SNR signals can be enhanced using the method outlined above, without 
unduly compromising the vertical resolution of the moving tool. It should be apparent that 
various combinations of prior art measurement methods and techniques can be used along 
with the time-domain averaging approach of the present invention. It should also be clear to 
those skilled in the art that the approach used in accordance with the present invention can 

15 also be used in measuring properties of samples of porous materials in a laboratory setting, 
as well as in situ logging-type including logging/measuring while drilling (LWD/MWD) 
measurements. The use of the method alone or in combination with other techniques, as 
discussed above is believed to be a significant contribution to the art of logging with wide 
ranging applications involving virtually all tools, and a broad range of practical 

20 applications, including logging while drilling. 

While the invention has been described with reference to the preferred 
embodiments, it will be appreciated by those of ordinary skill in the art that modifications 
can be made to the structure and form of the invention without departing from its spirit and 
25 scope which is defined in the following claims. 



30 



35 



-8- 



NY2- 1039344.1 



